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ABSTRACT 

Gaseous H2O has been detected in several cold astrophysical environments, where the observed 
abundances cannot be explained by thermal desorption of H2O ice or by H2O gas phase forma- 
tion. These observations hence suggest an efficient non-thermal ice desorption mechanism. Here, 
we present experimentally determined UV photodesorption yields of H2O and D 2 ice and deduce 
their photodesorption mechanism. The ice photodesorption is studied under ultra high vacuum con- 
ditions and at astrochemically relevant temperatures (18-100 K) using a hydrogen discharge lamp 
(7-10.5 eV), which simulates the interstellar UV field. The ice desorption during irradiation is mon- 
itored using reflection absorption infrared spectroscopy of the ice and simultaneous mass spectrom- 
etry of the desorbed species. The photodesorption yield per incident photon, Y pd (T, x), is identical 
for H2O and D2O and its dependence on ice thickness and temperature is described empirically by 
Ypd(T, x) — Y p d(T,x > 8)(1 — e~ x / l ( T )) where x is the ice thickness in monolayers (ML) and l(T) a 
temperature dependent ice diffusion parameter that varies between ~1.3 ML at 30 K and 3.0 ML at 
100 K. For thick ices the yield is linearly dependent on temperature due to increased diffusion of ice 
species such that Y pd (T, x > 8) = 10~ 3 (1.3 + 0.032 x T) UV photon" 1 , with a 60% uncertainty for 
the absolute yield. The increased diffusion also results in an increasing H2 0:OH desorption product 
ratio with temperature from 0.7:1.0 at 20 K to 2.0:1.2 at 100 K. The yield does not depend on the 
substrate, the UV photon flux or the UV fluence. The yield is also independent on the initial ice 
structure since UV photons efficiently amorphize H2O ice. The results are consistent with theoretical 
predictions of H 2 photodesorption at low temperatures and partly in agreement with a previous ex- 
perimental study. Applying the experimentally determined yield to a Herbig Ae/Be star+disk model 
provides an estimate of the amount of gas phase H 2 that may be observed by e.g. Herschel in an 
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example astrophysical environment. The model shows that UV photodesorption of ices increases the 
H2O content by orders of magnitude in the disk surface region compared to models where non-thermal 
desorption is ignored. 

Subject headings: Astrochemistry, Molecular data, Molecular processes, Methods: laboratory, Ultra- 
violet: ISM, ISM: molecules, Circumstellar matter 

1. INTRODUCTION 

H 2 0, in solid or gaseous form, is one of the most common species in molecular clouds, typically only second to H 2 
and sometimes to CO. This makes H2O, tog ether with CO, the dominant reservoir of oxygen during the critical stages 



of star formation 



van Dishoeck et al 



19931 ). H 2 is thus a key molecule in astrochemical models and its partitioning 



between the grain an d gas phase therefore has a large impact on the possible chemical pathways, including the formation 



Charnlev et al. 


1992: 


van Dishoeck 


2006) 



2008; 



Ioppolo et al 



2O forms througl 


1 hydr 


Tielens & Hagen 


1982 



hydrogenation of O (O2 or O3) on cold (sub)micron-sized silicate grain 



Leger et al 



1985 



Boogert fe Ehrenfreunc 



2004; 



Mivauchi et al 



2008). Other ices, like NH3 and CH4, probably form similarly, but observations show that H2O is 



of hydrogen nuclei. G as phase H2O formation is only efficient above 300 K (jElitzur fe Watson 



1978; 



Charnlev 



10 7 ( Bergin et al 



1 x 10 4 with respect to the number density 


( Elitzur & Watson 


1978 


Elitzur & de Jons 



1997 ). A t lower temperatures gas phase ion- molecule reactions maintain a low H2O abundance around 



1995). Any higher abundan ces require either thermal desorption of the H2O ice above ~100 K or 



non-thermal desorption at lower temperatures (jBergin et al 



1995 



Eraser et al 



20011) 



Gas phase H2O is observed from the ground only with great difficulty. Still both isotopi c and normal H2O have 



been detected in astrophysical environments from ground based telescop es (|Jacq et al 



Cernicharo et al 



1990; 



Gensheimer et al 



1996; 



1988; 



Knacke fe Larson 



1991 



van der Tak et al 



2006). The Infrared Spa ce Observatory (ISO) de 



tected warm H 2 Q gas unambiguously toward several low- and high-mass young s tellar objects ([van Dishoeck fe Helmich 



1996; 



Ceccarelli et al 



1999; 



Nisini et al 



1999 



Boonman fc van Dishoeck 



20031 ). ISO was followed by two other space 



based telescopes, the Submillimeter Wave Astronomy Satellite (SWAS) and Odin. In difference to ISO, SWAS and 



Melnick et al. 


2000; 


Hialmarson et al. 


2003) 



2003). Both telescopes have observed H2O gas toward star forming 



regions and detected the ex pected enhancements near protostars and in outflows where thermal ice desorption or gas 



phase formation is possible (jHialmarson et al 



2003 



Franklin et al 



also been detected toward photon dominated regions 



Snell et al 



2000; 



2008) . Critical for the pr esent study, H2O gas has 



Wilson et al 



20031 ) and is also more abun- 
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dant toward diffuse than toward dense clouds. These two results point to an efficient ice photodesorption mechanism 



(jMelnick fc Bergin 



2005) . The importance of photodesorption at the edges of clouds has more recently been modeled 
by Hollenbach et al. (ApJ, in press). They find that the H 2 gas abundance is enhanced by orders of magnitude at 
A y=2-8 mag into the clo ud when including photodesorption of H2O ice in their model at a rate derived from the results 



by 



be be large. 



Westlev et alj (Il995af). Circumste llar disks is a second region where the impact of photodesorption is expected to 
Willacy fc Langerl (|2000h showed that a photodesorption yield of 10 -3 molecules per photon is enough to 



explain observed gas-phase CO abu ndances in the outer r egio ns of flared disk s. Employing a similar photodesorption 
yield for H2O ice, the disk models of \ 



Dominik et al 



Willacvl (2007) both predict the existence of significant 



(120051) and] 

amounts of gas phase H2O in a layer above the midplane region. 
With the advent of the Herschel Space Observatory, co ld and warm H2O gas ob servations on scales of protostellar 



envelopes and disks will for the first time become possible 



van Kempen et al 



2008) . In preparation for these and other 



observations, and to interpret data from Odin and SWAS, non-thermal processes need to be better understood. These 
non-thermal desorption processes include ion/electron sputtering, desorption due to the release of chemical energy and 
photodesorption. Of these, sputt ering of ice by electrons and ions, ha s been investigated over a range of conditions 



during the past few decades (e.g. 



Brown et al 



1978 



Fama et al 



20081 ) and the dependencies of the yield on e.g. ice 



temperature, projectile type and energy are rather well understood. In contrast, only a handful o f laboratory studies 



e xists on the efficiency of ice photodesorption (e.g. 



Westlev et al 



Westlev et al 



1995a; 



Oberg et al 



2007b 



2008). 



(jl995allbl ) determined the photodesorption rate of H 2 ice experimentally to be 3 — 8 x 10 -3 molecules 
per UV photon for a 500 nm thick H2O ice. In their experiment the photodesorption rate depends on UV fluence 
as well as temperature. The photon fluence dependence, together with the observed gas phase H2 and O2 during 
irradiation, was taken as evidence that H2O photodesorption at low temperatures mainly occurs through desorption of 
photoproducts rather than of H2O itself. Several questions remain regarding the applicability of their study to astro- 
physical regions due to the uncertainty of the proposed mechanism. In addition, their dependence on photon fluence 
is not reproduced in recent CO and CO2 photodesorption experiments and cannot easily be e xplained theoretically 



(Andersson et al 



2006 



Oberg et al 



In a different experiment 



2007b; 



Yabushita et al 



Andersson fc van Dishoeck 



2008; 



Oberg et al 



20081) 



(2006) investigated H-atom photodesorption from H2O ice during irradia- 
tion at 157 and 193 nm using time-of-flight mass spectrometry. They found that the temperature and hence the origin 
of the desorbed H atoms varies significantly between crystalline and amorphous ice at 100 K. Thi s indica tes that pho- 



todesorption depends on the ice morphology, which is in contrast to the findings of 



Westlev et al 



(1995 



aj). Desorption 



2000). Both exper- 



of recombined D2 during irradiation of D2O ice at 12 K has also been found by (jWatanabe et al. 
iments provide valuable input for models, but cannot directly be used to determine the total H2O photodesorption 
rate. 
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UV irradiation of H2O ice results in photochemistry products as well as photodesorpt ion. A variety of photo- 



chemistry products are r eadily produced during irradiation of H2O dominated ice mixtures (jd'Hendecourt et al 



Allamandola et al 



Gerakines et al 



1996; 



1982 



1988 ). Pure H 2 Q ice pho tolysis results in the production of H2O2, OH ra dicals and HQ 2 a t 10 K 



Westlev et al 



Gerakines et al 



dl996h 



1995bt ). After a fluence of - 5 x 10 18 UV photons cm" 2 
found that the final band area of the formed H2O2 was only ~0.25% compared to the H2O band area and the OH 
band area was even smaller. No study exists fo r higher temperatures , but as photodissociation fragments become more 



mobile, e.g. O2 formation would be expected (jWestlev et al 



mm 



Only a handful mode l s of ice photodesorption exists in the literature. 



Andersson et al 



(2006); 



Andersson fc van Dishoeckl (|2008l ) have investigated H2O photochemistry and photodesorption theoretically using 
classical dynamics calculations. In the simulations they followed H2O dissociation fragments, after the absorption of 
a UV photon, in the top six monolayers of both crystalline and amorphous H2O ice at 10 K. For each ice they found 
that desorption of H2O has a low probability (less than 0.5% yield per absorbed UV photon) for both types of ice. 
The total H2O photodesorption yield from the top six ice layers was calculated to be ~ 4 x 10~ 4 molecules per incident 
UV photon. 

In the present study we aim at determining experimentally the photodesorption yields of H2O and D2O and their 
dependencies on ice thickness, temperature, morphology, UV flux and fluence as well as irradiation time. We use these 
results as input for an astrophysical model of a typical circumstellar disk to estimate the impact of photodesorption 
and to predict the observable column densities of H 2 as relevant to e.g. upcoming Herschel programs. 

2. EXPERIMENTS AND DATA ANALYSIS 

2.1. Experiments 



Fuchs et al. 


(2006 


) and 


Oberg et al. 


(2007b) 



The experimental set-up (CRYOPAD) is described in detail bv lFuchs et al. 
set-up allows simultaneous detection of molecules in the gas phase by quadrupole mass spectrometry (QMS) and in the 
ice by reflection absorption infrared spectroscopy (RAIRS) using a Fourier transform infrared (FTIR) spectrometer. 
The FTIR covers 1200 - 4000 cm -1 with a typica l spectral resolution of 1 cm -1 . The experimental procedure to derive 



photodesorption yields is described extensively in 



Oberg et al 



(2008), where the photodesorption measurements of CO, 



N2 and CO2 ices are reported. Here the procedure is summarized and only modifications to the procedure are described 
in detail. 

In the experiments, H2O and D2O ices of 1.5-28 monolayers (ML) are grown under ultra high vacuum conditions 
(P ~ 10 -9 mbar with the background pressure dominated by H2) at 18 - 100 K on a gold substrate that is mounted 
on the coldfinger of a He cryostat. The H2O sample is prepared from de- ionized H2O that is purified through at least 
three freeze-thaw cycles. The D2O sample is measured to have a 90% isotopic purity and is similarly freeze-thawed 
before use. 
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TABLE 1 

Summary of H2O and D2O experiments 



Experiment Composition Temperature [K] Thickness [ML] UV Flux [10 13 cm 2 s x ] 



1 


H 2 


18 


14 


2.3 


2 


H 2 


18 


15 


3.5 


3 


H 2 


100 


12 


1.1 


4 


H 2 


100 


13 


5.0 


5 


H 2 


100 


17 


2.3 


6 


D 2 


18 


10 


2.3 


7 


D 2 


18 


11 


3.5 


8 


D 2 


18 


17 


2.3 


9 


D 2 


30 


2.1 


3.5 


10 


D 2 


30 


3.2 


3.5 


11 


D 2 


30 


8.9 


3.5 


12 


D 2 


30 


11 


3.5 


13 


D2O 


40 


6 


1.7 


14 


D 2 


40 


14 


1.1 


15 


D2O 


40 


14 


3.5 


16 


D2O 


60 


9.6 


1.1 


17 


D 2 


60 


17 


3.5 


18 


D2O 


100 


1.5 


1.9 


19 


D 2 


100 


2.3 


1.8 


20 


D 2 


100 


5.1 


2.3 


21 


D2O 


100 


5.3 


2.3 


22 


D 2 


100 


6.8 


1.8 


23 


D 2 


100 


12 


5.0 


24 


D 2 


100 


13 


1.1 


25 


D2O 


100 


14 


3.5 


26 


D 2 


100 


16 


2.3 


27 


D 2 


100 


28 


1.7 


28 


D 2 O a 


30 


14 


2.3 


29 


H 2 8 


20 


14 


5.0 


30 


H 2 8 


30 


14 


5.0 


31 


H 2 8 


100 


14 


5.0 


32 


D 2 0/H 2 O b 


30 


9.5/48 


2.3 


33 


N 2 /D 2 O b 


18 


20/14 


2.3 



a Annealed at 100 K for 1 hour 
b Ice layers 
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Within the experimental uncertainties, we find that there is no difference in the photodesorption rate of 9.5 ML D2O 
ice deposited on top of 48 ML H2O ice compared with 8.9 ML D2O ice deposited directly onto the gold substrate. 
Since the nature of the substrate has no influence on the photodesorption, all other experiments are carried out with 
H2O or D2O ices deposited directly on the gold substrate. 

The ice films are irradiated at normal or 45° incidence with UV light from a broa dband hydrogen microwave dis charge 



2003). The 



lamp, which peaks around Ly a at 121 nm and covers 115-170 nm or 7-10.5 eV (jMunoz Caro fc Schuttd i 
lamp UV photon flux is varied between 1.1 and 5.0 xlO 13 photons cm~ 2 s _1 in th e different experime nts. The lamp 
flux is measured by a NIST calibrated silicon diode and actiometry as described in 



Oberg et al 



(2008) 



Table [T] summarizes the experiments in this study. The majority of the photodesorption experiments is carried out 
with D2O ice. These experiments are complimented with H2O and H3, 8 experiments to test isotope effects and to 
ensure the validity of mass spectrometric detections of OH (OD) and H2O (D2O) fragments. Layered experiments 
with H2O and D2O at 30 and 100 K and with N2 on top of D2O at 18 K are performed to check for substrate effects 
and to determine the ice loss behaviour when desorption is hindered. 

2.2. Data Analysis 

The UV induced ice loss rate during each H2O and D2O experiment is determined by RAIRS of the ice as a function 
of UV fluence. The intensity of the RAIRS profile is linearly correlated with the ice layer thickness up to ~20 ML, 
but the RAI RS profile can be u sed up to 50 ML for analysis as long as the non-linear growth above 20 ML is taken 



into account (jQberg et al 



2008). One monolayer is generally taken to consist of ~10 15 molecules cm' 2 and the loss 
yield, in molecules photon -1 , of the original ice is subsequently derived from the intensity loss in the RAIR spectra as 
a function of fluence. 

The determined ice loss yield is not necessarily the photodesorption yield. H2O has only dissociative transitions 
in the wavelength region of the lamp. Hence , the UV irradiation induces photodesorption as well as photochemistry 



( Gerakines et al 



1996; 



West lev et al 



1995a). UV irradiation may also induce structural changes in the ice that 
modify the infrared spectral features. These bulk processes, photolysis and rearrangement, are separated from the 
photodesorption by exploiting the different kinetic order behavior of bulk processes and surface desorption, i.e. first 
versus zeroth order processes. This method is described in detail in Oberg et al. (2008). 

Using RAIRS to determine the desorption of molecules depends on a reliable conversion between the ice infrared 
absorbance and the amount of ice molecules. Due to the fact that all ice measurements are do ne using RAIRS, th e 



Oberg et al 



( 200811 



ice thickness cannot be estimated from previously determined ice transmission band strengths. In 
the CO and CO2 appropriate RAIRS band strengths are reported. The H2O and D2O band strengths are estimated 
by assuming that the relative band strengths of CO, CO2 and H2O ice are th e same in tran s missio n and reflection- 
absorption spectra. This is found to be accurate within a factor of two by 



Ioppolo et al 



(2008). The thickness 
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uncertainty is then ~50%. For th e conversion between H2O and CO and CO2 band strength s , the measured band 



strengths of 



Hud gins et al 



(|1993[ ) are used after modification as suggested by 



band strengths of H 2 and D 2 were measured by 



Venvaminov k. Prendergast 



Boo gert et al 



([19971 ). The relative 



(|1997h . This results in H 2 and D 2 



stretching band strengths of 0.95 and 0.68 cm -1 ML -1 , respectively, for our set-up. These band strengths are converted 
to cm molecule -1 assuming a monolayer density of 10 15 molecules ML -1 cm -2 . 

Kinetic modeling of the integrated RAIRS profiles as a function of UV fluence and the determined band strengths 
together provide a total ice photodesorption yield. The simultaneous mass spectrometry of gas phase molecules during 
irradiation reveals the nature of the desorbed species, i.e. what proportion of H 2 ice desorbs as H 2 molecules versus 
photo-produced radicals and molecules. This is limited by the fact that less volatile molecules adsorb onto the heating 
shield and other semi-cold surfaces inside the experiment before reaching the mass spectrometer. Hence, the relative 
abundance of species with very different cryopumping rates, like H2 and H2O, cannot be derived. It is however possible 
to estimate the ratio of the predicted main desorption species: H2O and OH. 

The main sources of uncertainty in these experiments are the photon flux at the sample surface and ice thickness 
calibrations of ~30% and 50%, respectively. In addition, from repeated experiments, the H2O experimental results 
vary with 20%. The total uncertainty is hence ~60% for the total photodesorption rate. The relative desorption yields 
of different desorption products is more uncertain due to the additional assumptions that go into their derivation, i.e. 
that the products have similar pumping rates and QMS detection efficiencies. We estimate that the relative desorption 
yields thus have ^30% uncertainty in addition to the uncertainty of the total photodesorption yield. 

3. RESULTS 

3.1. Photodesorption process and products 

Figure Q] shows the RAIR spectrum of the D2O stretching band at 40 K as a function of UV fluence. At 40 K the 
UV photons simultaneously induce dissociation of bulk D2O ice and desorption of surface molecules. This is clearly 
seen in Fig. [51 where the ice thickness is plotted versus UV fluence. The ice loss is modelled by a combination of an 
exponential function and a lin ear function, where the exponential function describes the photolysis of bulk D2O as a 



first order process similarly to 



Cottin et al 



(2003). The linea r part of the ice los s is interpreted as photodesorption 



2007a)- To test that this model holds, a 



of surface molecules, which should be a zeroth-order process (jOberg et al. 
D2O ice is irradiated at 18 K when covered by a thick N 2 ice, which hinders desorption. The ice loss curve is then 
very well fitted (x 2 =l-8 for 13 data points) by the exponential function derived from a bare 18 K ice. This shows that 
the observed exponential decay is indeed due to bulk photolysis with a possible contribution from ice rearrangement. 
As discussed below the ice structure is affected by the UV irradiation, but this rearrangement seems complete within 
an UV fluence of 1 x 10 17 photons cm -2 . It may hence contribute somewhat to the exponential decay of the RAIRS 
profile, but is completely filtered out from the photodesorption rate determination. 



Oberg et al. 
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Fig. 1. — Spectra of the D2O stretching mode at 40 K as a function of UV fluence. The decreasing ice loss yield with fluence is due to 
the fact that ice is lost through a combination of bulk photolysis and photodesorption. The small bump at 2730 cm -1 is probably caused 
by the free OD stretch. 



There are no other species than H2O (D2O) visible in the RAIR spectra at any temperature. OH (OD) formation 
cannot be excluded, however, due to the spectral overlap of OH (OD) and H2O (D2O) transitions. Despite this overlap, 
photolysis of H2O into OH will result in a measured decrease of the H2O stretching band area. This is both because of 
the expected lower band strength of OH compared to H2O and because the rem aining H2O stretching band strength 



decre ases when the H2O network is disturbed by other molecules or fragments (| Oberg et al 



20071 ). The lack of H2O2 formation is in contrast to e.g. 



Gerakines et al 



2007a : 



Bouwman et al 



(1996). This is not a contradictory result 



however, since less than 0.1% of the H 2 ice is expected to be c onverted after a similar fluence, which is close to the 



1959). 



detection limit here for the strongest H2O2 band at 2850 cm -1 (jGiguere fc Harvevl l 

During irradiation H2 (D2) is always detected by the mass spectrometer. At the highest fluxes OH (OD) and H2O 
(D2O) are detected as well. A 14 ML H^O ice is irradiated at 20, 30 and 100 K while acquiring mass spectra (Fig. 
[3]) to quantify the relative desorption amounts of OH and H2O without any overlap with background H2O. The figure 
shows that the desorbing fraction of OH and H2O changes somewhat with temperature from 1.0:0.7 at 20 K to 1.2:1.4 
at 30 K to 1.2:2.0 at 100 K. No other species are observed at 20 and 30 K. In contrast, at 100 K, O2 is photodesorbed 
as well. It is important to note that O2 and H2O have very different cryopumping rates and hence their relative 
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Fig. 2. — The D2O ice layer thickness as a function of UV fiuence together with the fit (solid line), which is split up in ice loss due to 
bulk photolysis (dotted line) and photodcsorption (dashed line). The fitted bulk photolysis contribution is offset for visibility. 



mass spectromeric signals are not representative of their relative desorption rates. The upper limit of 2 desorption is 
estimated from the fact that a factor of 1.9 more OH and H2O is detected at 100 K compared to at 20 K, while the total 
photodcsorption rate from RAIRS increases by a factor of 2.4. Hence at most one fifth of the H 2 photodesorbs as 2 
at high temperatures. Desorption of H 2 2 cannot be excluded, even though it is not detected, since it is notoriously 
difficult to detect with a QMS. No build-up of H 2 2 is observed in the ice, which makes it unlikely that a large part 
of the ice is desorbed in the form of H 2 2 . Disregarding the small amount of H 2 that does not desorb as either OH 
or H 2 0, the H 2 yield relative to the total photodesorption yield is fitted linearly as a function of temperature for 
the 14 ML thick ice . This empirical fit yields an expression for the H 2 yield, Y pd ^ 2Ch as a function of the total 
photodesorption yield, F p( j: 



^pd,H 2 = /h 2 X Ypd (1) 

/h 2 o = (0.42 ± 0.07) + (0.002 ± 0.001) x T (2) 

/h 2 o + /oh ~ 1 (3) 

where / x is the fraction of the total photodesorption that occurs through species x. The relative yields are probably 



10 
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Fig. 3. — The mass spectra recorded during irradiation of a 14 ML thick H^O ice at 20, 30 and 100 K. In addition to photodesorbed 
ice product signals, background signals from H2O, OH, O and O2 are also seen. 



somewhat thickness dependent (jAndersson et ajj |2006) , but due to experimental constraints it is not possible to probe 
the relative yields for thinner ices. 

3.2. Yield dependencies on temperature, fluence, ice thickness, flux and isotope 
3.2.1. Temperature and photon fluence 

Figure [4] shows the combined bulk photolysis and photodcsorption for ^10 ML thick D2O ices between 18 and 100 K. 
The ice loss is dominated by bulk photolysis at low temperatures and by photodesorption at higher temperatures. This 
conclusion follows from the observed facts that 1) the degree of steady state photolysis decreases with temperature and 
2) that the photodesorption yield per incident photon, as measured by the slope of the linear part of the fit, increases 
with temperature. This increase in photodesorption yield with temperature is shown explicitly in Fig. [5^. Between 
18 and 100 K the dependence of the photodesorption yield on temperature for thick ices (>8 ML) is empirically fitted 
with a linear function: 



Yp d (T, x > 8) = 10~ 3 (1.3 (±0.4) + 0.032 (±0.008) x T) 
where T is the temperature in K and x the ice thickness in ML. The uncertainties are the model fit errors 



(4) 

- the total 
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Fig. 4. — Ice thickness (in ML) versus photon fluence (in 10 17 photons cm -2 ) for ~10 ML D2O ices at different temperatures, displaying 
the temperature dependent degree of ice bulk photolysis versus ice photodesorption. The curves are offset for visibility. At lower 
temperatures the bulk photolysis dominates the ice loss, while it is not visible above 60 K. Simultaneously the photodesorption rate, the 
slope of the linear part of the ice loss, increases with temperature. 



uncertainty of the yield is 60% as stated above. 
Figure 0] also shows that the onse t of photodeso r ption i s immediate, i.e. there is no fluence dependence, which 



is opposite to what was observed by 



West lev et al 



(|1995al ). This difference may be explained by a non- linear H2O 



freeze-out during the experiment, which is observed in this experiment (Fig. [6]) and probably present in the Westley 
experiment as well. Even under UHV conditions there is always some H2O present in the vacuum chamber. In our 
experiment this leads to an ice deposition rate of ~0.1 ML hr _1 at equilibrium, but up to 1 ML hr - 1 is deposited 



during the first hour after cool down, which is of the same order as the photodesorption rate presented in 



Westlev et al 



( 1995a)) . In their set-up it would not have been possible to separate this increased freeze-out rate from a lower desorption 
rate at the beginning of each experiment. This problem is circumvented here by using D2O for most experiments and 
by letting the H2O freeze-out reach equilibrium before starting the experiment. 

3.2.2. Substrate and ice structure effects 

In Fig. [5^ one of the two additional points at 30 K represents the yield from a 9.5 ML D2O ice deposited on top of a 
H 2 ice layer, showing that the substrate has no effect on the desorption yield. At 100 K this cannot be investigated 
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Fig. 5. — The D2O photodcsorption yield as a function of temperature (a), ice thickness (b), lamp flux (c) and isotope (d). In 
the temperature plot (a) the ices are ~10 ML thick. The two grey diamonds mark the desorption rate from a D20:H20 layered ice 
and an annealed ice. In panel (b) the thickness dependence is plotted and fitted for ices at 30 and 100 K. In panel (c) the ices are 
12—16 ML thick and irradiated at 100 K. The experiments marked with grey triangles in panel (d) are carried out with H2O instead of D2O. 



due to mixing of the two layered ices, but D2O experiments with different ice thicknesses indicate that the substrate 
has no effect at any temperature for ic es thicker than 2 ML (Fig. [5]d). 



The structure independence seen by 



Westlev et al 



(|1995al ) is also confirmed here for a 30 K i ce that is annealed at 



19811 ). and is subsequently 



100 K, until the spectra display a low frequency shoulder typical for crystalline ice (|Hagenl l 
cooled down (the second additional point at 30 K in Fig. [5^,). Figure [7] shows that this can be explained by the fact 
that the annealed ice returns to an amorphous state upon irradiation with less than 10 16 UV photons. At 100 K the 
irradiation does not yield amorphous ice, probably because the temperature is high enough for displaced molecules to 
diffuse back into a crystalline structure. 

3.2.3. Ice thickness 

At 100 K, where photodesorption completely dominates over bulk photolysis, there is no photodesorption yield 
dependence on ice thickness for D2O ices between 8 and 28 ML (Fig. [SJa) suggesting that the H2O ice photodesorption 
is only important in the top few layers. At 30 K the photodesorption yield is constant down to at least 3 ML. At both 
temperatures the thickness dependence is fitted by a function of the type 
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Fig. 6. — The build-up of H2O following cool down to 18 K (without UV irradiation), which is due to the small H2O contamination 
always present, also under ultra high vacuum conditions. The non-linear behavior is a result of the time required to reach steady-state 
between H2O freeze-out on the substrate, the desorption of adsorbed H2O from the chamber walls and the H2O pumping. 



Y(T, x) = Y pd (T, x > 8)(1 - e -*/'C T )) (5) 

where Y(T, x) is the thickness and temperature dependent photodesorption yield, Ypd(T, x > 8) the yield for thick ices 
at a certain temperature, x the ice thickness and l(T) is an ice diffusion parameter, whose origin is discussed further 
below. The IDL routine rap fit is used to fit the data with the results: Z(100 K) = 3.0 ± 1.0 ML and Z(30 K) < 2.7 ML, 
with the best fit of Z (30 K) = 1.3 ML shown in Fig. [5Jd. Thus photodesorption occurs deeper in the ice at 100 K 
compared to at 30 K by at least a factor of two. Extrapolating this to lower temperatures gives 



l(T) ~ 0.6 + 0.024 x T 



(6) 



3.2.4. Lamp flux 

The independence of the photodesorption yield on lamp flux found by 



Wes tlev et al 



(jl995al ) for higher fluxes 

(1 - 5 x 10 14 photons cm 2 s x ) is also seen in this study between 1.1 — 5.0 x 10 13 photons cm 2 s 1 for different 
temperatures and is shown for 100 K ices in Fig. [St. 
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Fig. 7. — Spectra of ices at 30 K (left) and 100 K (middle), and at 30 K, but deposited at 100 K (right), before onset of irradiation and 
after irradiation by ~10 16 and 10 17 photons cm -2 as indicated in the middle figure. Note the disappearance of the crystalline feature 
around 2460 cm -1 (marked with arrows) with increasing UV fluence for the annealed ice. 



3.2.5. H 2 versus D 2 

Within the experimental uncertainties, there is no difference between the H2O and D2O photodesorption yields at 18 
or 100 K (Fig. [5ji). It should be noted that at 18 K the photodesorption yield of H2O is highly uncertain, because the 
H2O freeze-out rate dominates over the photodesorption rate during the experiment. Nevertheless, these experiments 
support the direct applicability of above results, on D2O ices, to H2O ices. 

4. DISCUSSION 

4.1. H2O photodesorption mechanism 

The UV photodesorption mechanism of H2O ice does not depend on flux or substrate and is hence most likely due 
to direct absorption of UV light by H2O molecules, resulting in dissociation of the molecule into fragments with excess 
energy. The experiments show that once dissociated one of four different outcomes ensues: the dissociated fragments 
(i) photodesorb directly, (ii) recombine and photodesorb or kick out a surface H2O mol ecule, (hi) freeze out in the ice 



or (iv) recombine and freeze out in the ice. These are the same outcomes described in 



Andersson et al 



( 20061 ) when 



performing molecular dynamics simulations on H2O photodesorption. At high temperatures (100 K) the dissociated 
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OH or O fragments are also mobile enough to recombine to e.g. O2. 

The importance of the different pathways depends both on the temperature and where in the ice the molecule is 
dissociated. With increasing temperature the mobility of molecules, atoms and radicals is expected to increase. This, 
among other things, increases the reaction probability of OH and H. Assuming the same dissociation yield at all 
temperatures, the amount of radicals in the ice should then decrease with temperature since the increased mobility 
with temperature increases the recombination rate. This agrees well with the observed decrease in the steady state 
photolysis yield between 18 to 100 K (Fig. 2]). The increased recombination rate also explains the increasing H2O/OH 
gas phase ratio with temperature during photodesorption. Finally, the increased mobility of OH may also account 
for the observed photodesorbed O2 at 100 K and its absence at temperatures below 30 K. Although it cannot be 
excluded that the O2 forms also at lower temperatures and is thermally desorbed following formation above 30 K 



( Acharvva et al 



20071) . 

The increased mobility with temperature is also reflected in the ice thickness experiments, where the photodesorption 
occurs down to greater depths in the ice at 100 K compared to 30 K. The factor of two or larger penetration depth 
into the ice at 100 K compared to 30 K (Fig. [SJd) agrees well with the increase of the photodesorption yield between 
30 and 100 K (Fig. [5^). The increased photodesorption yield with temper ature is then most like ly due to an increased 



mobility rather than the overcoming o f react ion barriers, as suggested by 



The simulations of 



Andersson et al 



Westlev et al 



(|l995bh . 



(|2006D for an ice at 10 K indicate that photodesorption is only efficient in the 
top 2-3 layers for cold ices. At larger depths freeze-out of the dissociation products completely dominates. This is in 
excellent agreement with the results of the thickness dependent experiments at 30 K, where photodesorption is only 
important in the top 3 ML. The simulation is run at ps time scales, while the experiments cover several hours. The 
agreement for low temperatures between theory and experiment hence indicates that only short time scale processes 
matter for determining the photodesorption yield at temperatures below 30 K. At higher temperatures longer time 
scale processes, like thermal diffusion and desorption, increase in importance. 

As mentioned above, the depth at which photodesorption occurs increases with temperature, but still there is a 
certain ice depth where freeze-out of the recombined H2O is the only outcome. At 100 K the measurements are 
accurate enough to confirm that the mobility of the molecules following photodissociation and recombination is well 
described by a mean-free-path type model c x (1 — e~ x /'( T )), where x is the ice thickness, I the diffusion mean-free-path 
and c the maximum desorption yield for infinitely thick ice. This is also the case for CO2 photodesorption fragments 



( Oberg et al 



2008f ) and may hence be a universal feature for molecules that photodesorb following dissociation. This 
mean- free-path behaviour is best explained with the desorption of the energetic photodissocation product or recombined 
molecule itself, but does not exclude the other outcome of the molecular simulations, i.e. that H2O molecules also 



desorb indirectly from a kick of a H atom, which originates from photodissociation of another H2O molecule. 



16 Oberg et al. 

4.2. Comparison with previous experiments 

The results here agree on several important points with those of IWestley et al 



( 1995 J ib 



b). The maximum total 



photodesorption yields (H2O + OH) are the same within the reported uncertainties of both experiments. The de- 
termined photodesorption yield is likely a robust value, which is not significantly affected by different lamp spectral 
energy distributions, order of magnitude different UV fluxes and ice thicknesses. The present study also agrees with 



Westlev et al 



(|1995aU bl) on the photodesorption yield dependence on temperature and on the identification of the main 
desorption products - H2, O2 and H2O - with the one exception that we also detect OH. 

The apparent fluence dependence in the Westley experiment can be explained with H 2 freeze-out during the early 
stages of the experiment, especially since they mention a large H2O background pressure in their experiment. This 
is also in agreement with the mass spectrometer measurements of desorbed species show n in their paper (H2 a nd 



Westlev et al 



(1995a) to 



O2), which do not show any fluence dependent yields. The apparent fluence dependence led 
suggest that at low temperatures desorption occurs through reactions between O and OH. They subsequently claimed 
that while low temperature photodesorption occurs through photochemistry, high temperature photodesorption is a 
direct process. From the experiments here it is more likely that both low and high temperature photodesorption 
processes are dominated by direct photodesorption, but at high temperatures there is some additional desorption due 
to photochemistry of OH and O fragments. 
Ion sputtering o f ices has been more thoroughly investigated than ice photodesorption and recent experiments by 



Fama et al 



(2008) on the temperature dependence of ion sputtering of H2O ice suggest that the desorption mechanism 
is comparable for photodesorption and ion sputtering following the initial excitation by a photon or an ion. In 
particular, in both photodesorption and ion sputtering experiments, the desorption seems highly dependent on the 
formation and subsequent behavior of radicals and molecular products in the ice. More results are however required 
on e.g. the thickness dependence and the resulting desorption products during ion sputtering to make an actual 
comparison between the two processes. The absolute ion-sputtering yield of H2O depends on both the ion energy and 
ionic species, but it is generally a factor of 10 3 to 10 4 higher than photodesorption yields i.e. close to unity. When 
evaluating the importance of the two processes in an astrophysical setting, it is important to note that the ion flux in 
most regions is orders of magnitude lower than the UV flux. 

4.3. Astrophysical consequences 

This study shows that pure H2O ice photodesorbs directly or indirectly following fast intermediate photochemistry 
during which the the photodissociated fragments recombine. The mechanism hence does not depend on the photon 
flux level or on build-up of radicals in the ice. This means that the yield derived in the laboratory can be directly 
applied to astrophysical environments. Deep into clouds and disks the rate may be considerably reduced due to e.g. 
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CO ice covers. At the edges of clouds and disks, where other ices have not yet formed at large abundances, the rate 
for pure H2O ice is directl y applicable. Here , as a t e st case, it is applied to t he dis k surrounding a Herbig Ae/Be star 



using models developed by 



Dullemond et al 



(I2OO1I) 



Dullemond fc Dominikl (|2004 ) to fit the observed spectral energy 



distributions of these objects. In the model the physical disk model is static and the chemistry is kept very simple, 
only including H2O freeze-out, thermal and non-thermal desorption, and no gas phase chemical network except for the 
recombination of OH to form H 2 0. 

The model star has a mass of 2.5 M©, a radius of 2.0 R© and an effective temperature of 10500 K, typical of a 
Herbig Ae/Be star, and it emits a pure blackbody spectrum. The accompanying disk has a mass of 0.01 M Q , with 
an _R _1 surface density profile, and an outer radius of 300 AU. To avoid a sharp truncation, the surface density 
decreases as R~ 12 beyond the outer radius. The inner radius is set by a dust evaporation temperature of 1700 K. The 
radiation field and dust temp erature throughout the disk are calculated using the radiative transfer package RADMC 



( Dull emond fc Dominik 



2004) and the resulting disk is in vertical hydrostatic equilibrium, with a flaring shape. The 
gas temperature is set equal to the dust temperature. 

Gas phase H2O is initially distributed uniformly throughout the disk at a constant abundance of 1.8 x 10~ 4 x nn, 
where nu is the total number of hydrogen nuclei. This is somewhat artificial since H2O forms on grain surfaces, but 
if the model is run long enough (here to ~ 10 6 years) the final distribution will not depend on the initial distribution 
between the gas and grain. The H2O freezes out or adsorbs onto grain surfaces with the rate coefficient fc ac j s : 



fc ads = (4.55 x lO^WK-^V 1 ) n K Jj^ (7) 
where T g is the gas temperature and M the molecular weight of water. The numerical factor assumes unit stickin g 



efficiency, a mean grain radius of 0.1 /im and a grain abundance of 10 12 with respect to H2 (|Charnlev et al 
Once adsorbed onto the grains, the H 2 desorbs thermally with a rate coefficient k t hd- 



2001). 



«thd 



'1.26 x 10~ 21 cm 2 ) A—e 



(8) 



where n s is the number density of solid water and T gr is the grain temperature. The numerical factor assumes the same 
grain properties as in Eq. [7] and 10 6 binding sites per grain. The pre-exponentia l factor, A, and the binding energy, 



E^/k, are set to 1 x 10 30 cm 2 s 1 and 5773 K, respectively (jFraser et al 
rate coefficient fc p( j is defined as: 



2001). Finally the H 2 photodesorption 



fc pd = (3.14 x 10- 14 s- x ) G n H Y pd (9) 
where the numerical factor describes the UV photon flux onto a grain surface per unit time for the average interstellar 
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field (10 8 photons cm -2 s -1 ), Go is the scaling factor for the UV field that is output by RADMC for each grid point 
and Y p d is the photodesorption yield. I n addition to the external UV field a cosmic ray induced field is approximated 



by setting a lower limit on Go of 10 4 (jShen et al 



120041 ) . The photodesorption results in the release of both OH and 
H2O. The released OH is quickly rehydrogenated in the model in the gas, however, and hence we let all H2O ices 
desorb as H2O molecules with the yield: 



lpd,H a o(Tgr, <b) = 1(T 3 (1.3 + 0.032 x T gr ) (l - e -*/«D) (10) 

where x is the ice thickness in ML and l(T) the temperature-dependent diffusion length in ML (Eq. [JJ]). The model 
is run for two scenarios - (i) without and (ii) with photodesorption. Each run is 3 Myrs, the typical age of a Herbig 
disk, which results in steady state gas and grain phase abundances. Figure [8^ shows the gas phase H2O fraction 
in the disk without photodesorption as a function of radial and vertical distance from the central protostar. As 
expected the H2O is completely frozen out, except in the surface layer, when non-thermal desorption is excluded (Fig. 
Hk). When photodesorption is turned on, H2O is kept in the gas phase further in toward the mid-plane (Fig. |5Jd). 
Without photodesorption 0.6% of the H2O in the disk is in the gas phase at temperatures above the thermal desorption 
temperature of 100 K and 0.002% is in the gas phase below 100 K. With photodesorption included in the model, 0.6% 
of the H2O is still present as warm gas, but now 2% of the total H2O is in the gas phase at temperatures below 100 K. 
The total column density of warm T >100 K H2O gas, averaged over the entire disk, is hence the same in both cases 
(1.4 x 10 17 cm -2 ). In contrast the amount of cold H2O gas averaged over the entire disk increases from 5.0 x 10 14 to 
4.5 x 10 17 cm -2 when photodesorption is turned on. This means that a gas phase chemistry involving OH or H2O is 
possible deep toward the disk midplane also in the outer disk. For comparison we also run our model using a constant 
surface photodesorption yield of 10 -3 , which has often been used in the literature previously. For this particular model 
the total column density changes with less than a factor of two compared to using the derived yield equation from this 
study. The spatial distribution of gas phase H2O is however different using the different yields due to the fact that 
using the constant yield overestimates the desorption rate in the surface region and underestimates it deeper into the 
disk. 

This is a generic disk model commonly used to model disks around intermediate mass stars. To model an actually 
observed object would require a more detailed model that takes into account observed constraints on the disk structure. 
In addition, the calculated gas phase abundances may change somewhat when chemistry is taken into account. The 
general trend is however that photodesorption increases the amount of cold gas phase H2O by orders of magnitude. 
This is also the result of a recent PDR model showing that the inclusion of chemistry, while important for more 
accurate predictions, will not reduce the predicted column density of gas phase H2O dramatically (Hollenbach et al., 
ApJ, in press). Other effects such as grain growth may increase the photodesorption rate, but without a full chemical 
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Fig. 8. — Simulation of the distribution of gas phase over total H2O ratio in a circumstellar disk without (a) and with (b) photodesorption. 
The white, H^O-gas-dominated area extends more than 50 AU deeper into the disk when photodesorption is included, illustrating the 
large impact of photodesorption on the chemistry in the outer parts of disks. 

network it is unclear how much of this increase will be off-set by photodissociation of the desorbed H2O. Probably 
only observations of cold H2O on the scale of protoplanetary disks will yield an answer. The beam of the imminent 
Herschel Space Observatory is of the same order of magnitude as this modelled disk and hence these results show that 
large amounts of cold H2O will be observable. For exactly this purpose - to observe cold H2O gas - the WISH program 
was approved as a Herschel key program. The amount of gas phase H2O in disks due to non-thermal desorption may 
hence be answered very soon indeed. 



5. CONCLUSIONS 



20 Oberg et al. 

1. The total D2O and H2O photodesorption yields are indistinguishable within the experimental uncertainties 
and are empirically described by Ypd(T, x) = Y p d(T, x > 8)(1 — e ~ x / l (. T )^ w here Y p d(T, x) is the thickness 
and temperature dependent photodesorption yield, x the ice thickness in monolayers and 1{T) an ice diffusion 
parameter that varies between 1.3 ML at 30 K and 3.0 ML at 100 K. 

2. For thick ices (>8 ML), the yield depends linearly on temperature such that Y p a(T, x > 8) = 
10~ 3 (1.3 + 0.032 x T) photon" 1 . The yields agrees, within the reported 60% uncertainty, with a previous ex- 



periment (W est lev et al 



1995a). 



3. The nature of the desorbed species is temperature dependent, with equal amounts of OH and H2O detected at 
low temperatures. At higher temperatures the H^OiOH fraction is ~ 2:1 and in addition about a fifth of the ice 
photodesorbs as heavier fragments like O2. The fraction of the total photodesorption yield that results in H 2 
molecules desorbing is described by /h 2 o = (0.42 ± 0.07) + (0.002 ± 0.001) x T. 

4. We find no yield dependence on p hoton flux or fluence. The fluence independence is in contrast with a previous 



experiment (jWestlev et al 



mm - 



5. We also find no dependence on the ice structure i.e. whether the D2O ice is amorphous or crystalline. This is 
consistent with spectroscopic evidence of fast destruction of crystalline ice into an amorphous state following UV 
irradiation. 

6. The photodesorption yield and dependencies found here are consistent with previous theoreti cal predictions of 



H2O photodesorption, where the photodesorption is limited to the top few layers of the ice (jAndersson et al 



2006) . In addition we sec that the photodesorption yield increases with ice temperature because of the increased 



mobility of the photolysis fragments, allowing desorption from deeper within the ice. 

7. Applying the experimental yield to a Herbig Ae/Be star+disk model we calculate that the predicted amount 
of cold (< 100 K) gas phase H 2 0, averaged over the entire disk, increases with orders of magnitude due to 
photodesorption. 
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